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Absorption spectroscopy is routinely used to characterise chemical and biological samples. For 
the state-of-the-art in absorption spectroscopy, precision is theoretically limited by shot-noise 
due to the fundamental Poisson-distribution of photon number in laser radiation. In practice, 
the shot-noise limit can only be achieved when all other sources of noise are eliminated. Here, 
we use wavelength-correlated and tuneable photon pairs to demonstrate sub-shot-noise absorp¬ 
tion spectroscopy. We measure the absorption spectra of spectrally similar biological samples— 
oxyhaemoglobin and carboxyhaemoglobin—and show that obtaining a given precision in resolution 
requires fewer heralded single probe photons compared to using an ideal laser. 


A sample’s absorption spectrum is typically measured 
by comparing the wavelength and intensity of incident 
light with the wavelength and intensity of transmitted 
light. However, the state-of-the-art is bound in precision 
by the shot-noise limit (SNL) which can ultimately limit 
precision in practice due to acceptable optical effects on 
the sample itself from the probe, including damage Hl- 
13:. Here we demonstrate a statistical benefit to using 
frequency-correlated photon pairs when performing ab¬ 
sorption spectroscopy. By using heralded single photons, 
there is lower fundamental noise than that of ideal laser 
emission of equal intensity. This is advantageous in spec¬ 
troscopy where the error in measurement is optical-power 
dependent, such as in Doppler thermometry mm, and 
for measuring with great precision in short time inter¬ 
vals, while minimising photo-chemistry, such as for ob¬ 
serving cell dynamics [2J. We measure the spectra of two 
different types of blood protein— haemoglobin bound to 
oxygen and haemoglobin bound to carbon monoxide— 
with sub-SNL precision per detected photon. Resolving 
two different samples in this way demonstrates capacity 
to achieve a greater resolution of distinguishing two ab¬ 
sorption features in the same sample. We support this 
by analysing sub-SNL performance when measuring an 
optical filter. 

Correlated pairs of photons can be used to herald the 
generation of single-photons [0 which can in turn be used 
to measure optical transmission with precision beyond 
the SNL [7,:8] with optimal performance [91. The increase 
in precision over an attenuated laser comes from the re¬ 
duction of possible photon number detection outcomes 
for correlated photon pairs. Detecting one of the photons 
in each pair heralds the presence of its counterpart—a 
single probe photon—which is either absorbed or is not. 
By contrast, coherent and thermal sources of light con¬ 
tain more than vacuum and single photon contributions, 
even when attenuated —this leads, at best, to Poisson dis¬ 
tributions of the output photon-number detection statis¬ 
tics and error on any subsequent estimate. 

The performance of various probe states have been 


compared using estimation theory for various types of de¬ 
tection schemes Gang. Fock states | N) have been found 
to be optimal as probes since they have zero uncertainty 
of photon-number at the input, resulting in a minimum 
spread for the photon number at the output [7J £1 ■ In¬ 
creasing N decreases uncertainty in measurement, but 
does not provide additional advantage in scaling. This 
is quantitatively similar to the sub-SNL advantage that 
is obtainable in entanglement-enhanced phase estimation 
once losses and decoherence have been accounted for m 

Previous experiments that achieve sub-SNL precision 
to reduce photo-damage to samples mm equate to opti¬ 
cal phase estimation which relies on non-classical inter¬ 
ference via entanglement and squeezing. The quantum- 
enhancement pursued here is based on the sub-poissonian 
statistics of single photons [7: that can be achieved 
equally with single photon emitters or heralded pho¬ 
tons generated from parametric processes. This has 
been achieved for monochromatic absorption imaging 
using spatially correlated photon pairs [Sj. Exploiting 
frequency correlations, photon pairs have been used to 
reveal absorption spectra mHnj — for low illumination, 
correlated photon pairs enable a high signal-to-noise ra¬ 
tio to be maintained despite increasing levels of back¬ 
ground illumination Eg. However, these spectroscopy 
experiments did not achieve sub-SNL precision. With 
increased system efficiency, source brightness and exten¬ 
sion of spectral range, the results we present open the way 
to practical sub-SNL absorption spectroscopy for broad 
application across science and technology. 

Absorption of a sample, 0 < a < 1, is determined by 
comparing a known input intensity N and a measured 
output intensity N' = (1 — a)N. This is equivalent to es¬ 
timating the overall loss a = 1 — (1 — a 3 )(l — a 2 )(l—aq), 
where photon loss occurs (i) before the sample at the 
source (aq); (ii) due to absorption by the sample ( 0 : 2 ); 
(iii) after the sample during measurement (0:3). This 
treatment is valid, provided aq and a 3 are calibrated 
with high accuracy separately. For a dilute solution, «2 
is related to the absorbance, A, of the sample through 
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FIG. 1. Theory for precision in estimating a: use of 

the Fock state |1) is represented by the solid line and use of 
a coherent state with the same average intensity of IV = 1 
photons is represented by the dashed line. Note that diver¬ 
gence at a = 0,1 corresponds to estimates with vanishing 
variance—e.g. for a — 1, repeated trials with |1) will always 
yield a zero-photon detection event yielding A 2 a = 0. Inset: 
The dot-dashed line corresponds to the ideal quantum advan¬ 
tage that can be obtained using 11), defined by the precision 
ratio of A 2 a for the jl) and a coherent state with N = 1. 


the Beer-Lambert Law, which gives an exponential at¬ 
tenuation depending on the molar absorption coefficient, 
e, the molar concentration, c, and the distance through 
the sample, l: A = eel = ln[ 1 — 02 ] |19| . 

Precision in estimating a —the reciprocal of the root 
mean square error A a —is limited by the fundamental 
statistical fluctuations of the input probe due to the 
quantum nature of light [ 201 . Given v repetitions of 
measurements on the output state p' for the output pho¬ 
ton number, the precision for estimating a is given by 

An = A p 'N / y/vN, where A p >N = ^J ( N 2 ) p / — {N) 2 p/ , and 

N is the output photon-number operator. Absorption 
spectroscopy currently uses laser emission as its lowest 
noise input probe, and so the probability to detect n 
photons from the input probe, P(n), and n! photons after 
absorption, P'(n') : are both governed by a Poisson dis¬ 
tribution: P (n) 1 —>■ P' ( n') = e ~N(i~a) [_/V(l — a)] n /n'\. 
This yields the fundamental limit of precision for clas¬ 
sical light (the SNL) A a c i = \/(l — a)/vN. Other 
incoherent broadband probe sources such as sodium 
lamps, and non-ideal laser setups lead to noisier in¬ 
tensity measurements. For a Fock state \N) acting as 
the probe, the loss process results in a binomial distri¬ 
bution P' (n') = (^) (1 — a)' 1 a N ~ n which results in 
A ap = y/a(l — a)/uN , an improvement of 1/i/a: over 
the SNL. Equivalently, for a given target precision Aa, a 
factor of (1 — a) fewer single photons in the state | 1 ) are 
needed than when using a coherent state with N = 1. 

We compare in Fig. [T] the performance of loss esti¬ 
mation for using the Fock state |1) against using an 
ideal laser modelled with Poisson-distributed photon- 
number statistics and having the same photon intensity, 


N = 1. Performance is quantified using zz/A 2 a which 
corresponds to the statistical information gained per de¬ 
tected photon—this is the same as the Fisher Informa¬ 
tion that is widely used in quantum parameter estima¬ 
tion mm- f^. |T] illustrates that in principle there is an 
advantage for using |1) for any a, however the magnitude 
of the improvement scales with a itself: the greatest im¬ 
provement occurs at low total absorption. Since a is de¬ 
fined including losses throughout the system, this scheme 
can obtain a quantum advantage with non-perfect com¬ 
ponents. 

Our experimental setup is shown in Fig. |2ja). A con¬ 
tinuous wave diode laser at 403.9 nm pumps a 10 mm 
long periodically poled potassium titanyl phosphate 
(PPKTP) non-linear optical crystal, phase-matched for 
collinear type-II spontaneous parametric down conver¬ 
sion (SPDC), to generate orthogonally-polarised photons 
in the same spatial mode (H and V). These are then 
separated deterministically using a polarising beamsplit¬ 
ter, while the blue pump laser is removed using three 
dichroic mirrors normal to the beam, two dichroic mirrors 
at 45°, one longpass filter and one wide bandpass filter 
centred at 810nm (FWHM= 50 nm). The photon pairs 
are generated in a low gain parameter regime of SPDC 
and so are modelled by | SPDC) ~ |0) a |0)b + e|l) a |l)j,. 
The crystal is in a temperature-tuneable oven which we 
use to control the wavelength correlations of the emit¬ 
ted photons—this is pre-calibrated using a single-photon 
sensitive spectrometer (see inset of Fig. [ 2 ]). This con- 
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FIG. 2. Setup for sub-shot-noise spectroscopy, high¬ 
lighting the simplicity of the scheme. Wavelength cor¬ 
related photon pairs are generated using a laser to pump a 
non-linear optical crystal (PPKTP) that is phase-matched 
for collinear type-II SPDC and temperature tuned for wave¬ 
length control. The sample to be measured can include opti¬ 
cal filters, cuvettes containing a liquid chemical or biological 
sample. PBS-polarising beamsplitter; MMF-multimode fibre; 
APD-avalanche photodiode. Inset. The calibrated temper¬ 
ature dependent joint spectrum of the horizontally polarised 
(blue) and vertically polarised (purple) photons generated in 
our PPKTP crystal when pumped with a 403.9 nm continu¬ 
ous wave diode laser. 
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figuration removes the need to optically determine the 
wavelength of either photon at the output of our setup 
during a sample measurement and simplifies the setup to 
use only two single photon detecting pixels—in our case 
single photon counting modules and coincidence logic. 

Absorption estimates from each single probe photon 
are distributed across the photon’s bandwidth -this is a 
limiting factor for both the resolution and precision of 
measured spectra. Here, the bandwidth of photons gen¬ 
erated from SPDC is dependent upon that of the pump 
laser. We use a laser with a linewidth of 0.057 nm 
(FWHM) resulting in linewidths of ~ 0.5 nm and ~ 
0.7 nm for the photons in arms a and b respectively. The 
accessible wavelength range for each arm for our setup 
is 773 < X a < 809 and 806 < At < 845 nm, restricted 
by the maximum temperature of our oven (= 200° C) and 
our use of a fixed frequency pump laser. Each path a and 
b is coupled into multimode fibres to increase photon col¬ 
lection efficiency and to reduce sensitivity to mechanical 
vibrations. In the absence of a sample, the efficiency of 
each arm is r] a ~ 35% and 777 , ~ 29%—this includes the 
~ 69% specified efficiency of our single-photon detector 
modules that are used with coincidence logic to record 
both the singles count rates (N a , N b ) and the coincidence 
count rates ( N ab ). 

Fig-i displays the application of our setup to mea¬ 
sure the spectral response of a Gaussian bandpass fil¬ 
ter []] placed in path b. The absorption is estimated di¬ 
rectly [2T] from the ratio of heralded single photon de¬ 
tection events in path a to coincidences, according to 

Ctexp = 1 ^ N ab /N a . (1) 

To take a complete spectrum, a exp is calculated at a 
range of different crystal temperatures corresponding to 
a known set of wavelengths. For each probing wavelength 
we compute the mean value of a exp = 1 — (1 — oq)(l — 
tt 2 )(l — 013 ) and the variance A 1 2 a exp over v = 1500 trials 
using a 1 s integration time per trial. To ensure a sta¬ 
ble probing wavelength after each temperature change, 
we enforce a 300 s stabilisation period before beginning 
the next set of measurements. The absorption spectrum 
of the filter, 0 : 2 , was found by dividing out the system 
absorption due to oq and «3 that was characterised sep¬ 
arately. The mean values of a 2 for each wavelength are 
plotted in Fig. [3ja) where the black dots are from our 
experimental setup and the black line is that from a clas¬ 
sical scan using a UV/Vis spectrometer. Close agree¬ 
ment verifies comparable accuracy to commercial spec¬ 
troscopy, up to 0.7nm resolution defined by the photon 
bandwidth in path b. There was a uniform 0.375nm off¬ 
set between our setup and the spectrometer due to cali- 
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FIG. 3. Sub-shot-noise absorption spectra of a con¬ 
trol feature. A Gaussian bandpass filter was placed in 
arm b and its spectral response across the wavelength range 
807 — 818 nm was measured with correlated photon pairs 
(dots) and a commercial UV/Vis spectrometer (solid line), 
(a) The mean absorption of the filter extracted from 1500 
estimates for each measured wavelength, (b) Measured quan¬ 
tum advantage (data points) for estimating a of the Gaussian 
filter, quantified as a ratio of the precision obtainable with an 
ideal laser (solid line). The theoretical maximum that ac¬ 
counts for our system efficiencies, using ideal quantum states, 
is computed using a exp and represented by the dot-dashed 
line. Error bars computed as explained in the main text. 


bration discrepancy between the UV/Vis and the single¬ 
photon-sensitive spectrometers, and minor alignment er¬ 
ror of the Gaussian filter when mounted in the UV/Vis 
spectrometer. 

The quantum advantage obtainable with our setup to 
measure the Gaussian filter is quantified in Fig. !b). 
We plot the ratio of our computed variance for the her¬ 
alded single photons, compared to what we would obtain 
using a shot-noise limited ideal laser (solid line in the in¬ 
set of Fig. m), for each corresponding mean estimate of 
a exp in FigFWa). Of the 1500 trials we use sets of 100 to 
compute 15 quantum advantage parameters; the mean of 
these are plotted in Fig. [3])b) and the standard error of 
each mean value are plotted as the error bars. The blue 
region shows where a quantum advantage is obtained and 
the dot-dashed line shows the maximum theoretical ad¬ 
vantage for each a exp value. A quantum advantage in 
precision is achieved across the range of the filter, with a 
maximum advantage of 22.20 ± 0.04% per detected pho¬ 
ton. 

We demonstrate sub-SNL spectroscopy of biologi¬ 
cal samples in the near-infrared region by measur¬ 
ing the absorbance spectra for two different types of 
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FIG. 4. Sub-shot-noise absorbance spectra of Hb0 2 and HbCO. Samples of Hb0 2 (red) and HbCO (brown) were placed, 
in turn, in path a and the spectral response was measured over the wavelength range ~ 790—808 nm. The experimental quantum 
advantage for estimating absorption a is plotted separately for Hb0 2 (a) and HbCO (b). The error bars are computed using 
the same treatment as for the Gaussian filter, (c) The mean value of a was converted to absorbance, A. (d) Resolving otHbO 
from etHbeo at 790.5 nm probe wavelength. For each total number of photons, 800 estimates for a are made and the mean is 
plotted (circles). The error on these estimates is quantified by one standard deviation (Aa) of these estimates and plotted as 
error bars. Solid black vertical lines represent the total number of photons detected in our experiment and required to resolve 
one a from each other by 2A a, 3Aa and 4A«—these lines are computed using fits to the error bars for each sample. We 
compare this directly to the total number of photons required to resolve the two a when using an ideal laser (vertical dashed 
lines). The shaded regions (left to right) correspond to 1250, 2850 and 4810 fewer required photons. We performed a high 
resolution scan (1ms increments) for the region of 2Aa separation. 


haemoglobin: oxyhaemoglobin (Hb 02 ) and carboxy- 
haemoglobin (HbCO). The samples were placed in path 
a and subjected to the same intensity and data integra¬ 
tion time as when measuring the Gaussian filter. From 
the mean estimate for 0 : 2 , we used the Beer-Lambert law 
(A = ln[ 1 — £* 2 ]) to calculate the estimated absorbance 
spectrum for each sample, plotted in Fig |4](c) . A near- 
flat spectrum with carboxyhaemoglobin more absorbing 
than oxyhaemoglobin is obtained, as expected both from 
literature |22j 123] and from our own characterisation us¬ 
ing the UV/Vis spectrometer. We plot in Fig 0a, b) 
the quantum advantage in precision of estimating the 
spectral profile of each sample, showing that a quantum 
advantage in precision per detected photon is achieved 
across the entire spectral range for both samples. 

Fock states can be used to reduce the number of pho¬ 
tons required to discriminate between different absorp¬ 
tions a , leading to a higher absorptive resolution than 
when using an ideal laser. We demonstrate this in 
Fig 0d) using the two different haemoglobin samples. 
Using one probe wavelength (790.5 nm) and a constant 
intensity, we computed a for an increasing total num¬ 
ber of photons—we controlled this by increasing the in¬ 


tegration time linearly starting from 1 ms, using 1 ms 
increments for high resolution and 5ms increments for 
lower resolution. For each increment we computed 800 
estimates of a and calculated the mean and standard 
deviation Aa. The standard deviation decreases with in¬ 
creasing total number of photons at a rate that is faster 
than can be obtained with an ideal laser. This is shown 
by the vertical lines in Figj4] (d) that quantify the num¬ 
ber of photons each scheme requires to resolve ctHbO from 
&HbCO by a separation defined by multiples of the stan¬ 
dard deviation, for our experiment (solid) and an ideal 
laser (dashed). The difference in total number of photons 
required in our scheme compared to using an ideal laser 
will depend on the values of a (c.f. Fi# 

We have demonstrated using correlated photon pairs 
for sub-SNL precision in absorption and absorbance spec¬ 
troscopy. There are three main aspects of our setup 
where performance can be improved, (i) Increasing over¬ 
all system efficiency will enable a greater quantum advan¬ 
tage (Fig. 0— superconducting detectors are a promis¬ 
ing approach, > 90% efficiency have been demonstrated 
at telecommunication wavelengths [24] with promise of 
high efficiency at shorter wavelengths [25] ; high efficiency 
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single-photon-sensitive cameras [5] could be applied alle¬ 
viating the need for cryogenic temperatures, (ii) Extend¬ 
ing emission to a larger spectral range would increase 
application—this can be accomplished with a tuneable 
pump laser [25]. (in) Generating photons with narrow 
bandwidths would enable higher spectral resolution -for 
example, cavity-enhancement T] would enable applica¬ 
tion to atomic spectroscopy. 

Sub-SNL correlated photon pair spectroscopy is a 
promising alternative to current forms of spectroscopy— 
where in practice most systems do not achieve SNL 
precision—and those where intensity fluctuations have 
been suppressed [27ff29] , Likely applications include 
characterising samples that are photo-sensitive or of low 
concentration or volume—the principle demonstrated in 
Fig. gd) applied to neighbouring points on a single 
spectrum, demonstrates the advantage when probing 
shallow spectral features. Unlike other demonstrations 
in quantum-enhanced precision measurement, our setup 
does not require entanglement or multiphoton interfer¬ 
ence. Equivalent performance could be obtained with 
extremely narrow bandwidth single photons generated 
from quantum dots E3DJEII or atom-cavity systems [32j . 
Proposed characterisation of photon-counting and homo¬ 
dyne detectors with Fock states [33] could use an iteration 
of our setup to benchmark spectral response. 

Acknowledgements: The authors would like to 
thank R. Oulton, A. Young, P. Androvitsaneas, S. Car- 
swell, B. Lang and E. Harbord for use of the spectrom¬ 
eter; S. Bellamy, P. Dunton and M. Mercer for wet-lab 
support; G. Magro for materials-lab support; X.Q. Zhou, 
J. Rarity, M. Padgett, A. Datta, for useful discussions; 
N. Tyler, A. Santamato and C. Wilkes for preliminary 
work on the source. The authors are grateful for sup¬ 
port from DSTL, EPSRC, ERC, NSQI. JLOB acknowl¬ 
edges a Royal Academy of Engineering Chair in Emerging 
Technologies. JCFM is supported by a Leverhulme Trust 
Early Career Fellowship. 


* Beccie.Whittaker@Bristol.ac.uk 
t Jonathan.Matthews@Bristol.ac.uk 

[1] F. Wolfgramm, C. Vitelli, F. A. Beduini, N. Godbout, 
and M. W. Mitchell, Nature Photonics 7, 28 (2013). 

[2] M. A. Taylor, J. Janousek, V. Daria, J. Knittel, B. Hage, 

H. -A. Bachor, and W. P. Bowen, Nature Photon. 7, 229 
(2013). 

[3] M. A. Taylor and W. P. Bowen, arXiv: 1409.0950. 

[4] T. Stace and A. Luiten, Physical Review A 81, 033848 

( 2010 ). 

[5] T. M. Stace, G.-W. Truong, J. Anstie, E. F. May, and 
A. N. Luiten, Phys. Rev. A 86, 012506 (2012) 

[6] P. Grangier, G. Roger, and A. Aspect, Europhys. Lett. 

I, 173 (1986). 

[7] E. Jakeman and J. G. Rarity, Opt. Comm. 59, 219 


(1986). 

[8] G. Brida, M. Genovese, and I. Rou Berchera, Nature 
Photon. 4, 227 (2010). 

[9] G. Adesso, F. Dell’Anno, S. De Siena, F. Illuminati, and 

L. A. M. Souza, Phys. Rev. A 79, 040305R (2009). 

[10] M. M. Hayat, A. Joobeur, and B. E. A. Saleh, J. Opt. 
Soc. Am. A , 16 (1999). 

[11] V. D’Auria, C. de Lisio, A. Porzio, S. Solimeno, and 

M. G. A. Paris, J. Phys. B: At. Mol. Opt. Phys. 39, 1187 
(2006). 

[12] A. Monras and M. G. A. Paris, Phys. Rev. Lett. 98, 
160401 (2007). 

[13] B. Escher, R. de Matos Filho, and L. Davidovich, Nature 
Physics 7, 406 (2010). 

[14] G. Scarcelli, A. Valencia, S. Gompers, and Y. Sliih, Appl. 
Phys. Lett. 83, 5560 (2003). 

[15] A. A. Kalachev, D. A. Kalashnikov, A. A. Kalainkin, 
T. G. Mitrofanova, A. V. Shkalikov, and V. V. Samart- 
sev, Laser Phys. Lett. 5, 600 (2008). 

[16] O. Slattery, L. Ma, P. Kuo, Y. S. Kim, and X. Tang, 
Laser Phys. Lett. 10, 075201 (2013). 

[17] A. Yabushita and T. Kobayashi, Phys. Rev. A 69, 013806 
(2004). 

[18] A. A. Kalachev, D. A. Kalashnikov, A. A. Kalainkin, 
T. G. Mitrofanova, A. V. Shkalikov, and V. V. Samart- 
sev, Laser Phys. Lett. 4, 722 (2007). 

[19] F. Miller, A. Vandome, and J. McBrewster, Beer- 
Lambert Law (VDM Publishing, 2009). 

[201 V. Giovannetti, S. Lloyd, and L. Maccone, Science 306, 
1330 (2004). 

[21] G. Brida, M. Genovese, and C. Novero, J. Mod. Opt. 
47, 2099 (2000). 

[22] W. Zijlstra, A. Buursma, and W. Meeuwsen-Van der 
Roest, Clinical chemistry 37, 1633 (1991). 

[23] W. Zijlstra and A. Buursma, Comparative Biochemistry 
and Physiology Part B: Biochemistry and Molecular Bi¬ 
ology 118, 743 (1997). 

[24] F. Marsili, V. Verma, J. Stern, S. harrington, A. Lita, 
T. Gerrits, I. Vayshenker, B. Baek, M. Shaw, R. Mirin, 
and S. Nam, Nature Photonics 7, 210 (2013). 

[25] F. Marsili, F. Bellei, F. Najafi, A. Dane, E. Dauler, 
R. Molnar, and K. Berggren, Nano Letters 12, 4799 
( 2012 ). 

[26] R.-B. Jin, R. Shimizu, K. Wakui, H. Benichi, and 
M. Sasaki, Optics Express 21, 10659 (2013). 

[27] S. Kasapi, S. Lathi, and Y. Yamamoto, Journal of the 
Optical Society of America B 17, 275 (2000). 

[28] E. S. Polzik, J. Carri, and H. J. Kimble, Phys. Rev. Lett. 
68, 3020 (1992). 

[29] S. Kasapi, S. Lathi, and Y. Yamamoto, Opt. Lett. 22, 
478 (1997) 

[30] R. Young, R. Stevenson, A. Shields, P. Atkinson, 
K. Cooper, and D. Ritchie, Journal of Applied Physics 
101, 081711 (2007). 

[31] H. Jayakumar, A. Predojevic, T. Huber, T. Kauten, 
G. Solomon, and G. Weihs, Physical Review Letters 110, 
135505 (2013). 

[32] A. Kuhn, M. Hennrich, and G. Rempe, Phys. Rev. Lett. 
89, 067901 (2002). 

[33] M. Barbieri, A. Datta, T. J. Bartley, X.-M. Jin, W. S. 
Kolthammer, and I. A. Walmsley, arXiv:1502.00681. 




